arXiv:1502.01388vl [astro-ph.SR] 4 Feb 2015 


THE MILKY WAY PROJECT: WHAT ARE YELLOWBALLS? 


C. R. Kerton 

Department of Physics & Astronomy, Iowa State University, 12 Physics Hall, Ames, lA 50011, 

USA 

kertonOiastate.edu 
G. Wolf-Chase^ 

Astronomy Department, Adler Planetarium, 1300 S. Lake Shore Drive, Chicago, IL 60605, USA 

gwolfchase@adlerplanetarium.org 
K. Arvidsson 

Trull School of Science and Mathematics, Schreiner University, 2100 Memorial Blvd., Kerrville, 

TX 78028, USA 

KDArvidssonSschreiner.edu 
C. J. Lintott and R. J. Simpson 

Oxford Astrophysics, Denys Wilkinson Building, Keble Road, Oxford 0X1 3RH, UK 

cjISastro.ox.ac.uk, robert.simpson@astro.ox.ac.uk 

ABSTRACT 

Yellowballs are a collection of approximately 900 compact, infrared sources identified 
and named by volunteers participating in the Milky Way Project (MWP), a citizen- 
science project that uses GLIMPSE/MIPSGAL images from Spitzer to explore topics 
related to Galactic star formation. In this paper, through a combination of catalog cross¬ 
matching and infrared color analysis, we show that yellowballs are a mix of compact 
star-forming regions, including ultra-compact and compact H II regions, as well as 
analogous regions for less massive B-type stars. The resulting MWP yellowball catalog 
provides a useful complement to the Red MSX Source (RMS) survey. It similarly 
highlights regions of massive star formation, but the selection of objects purely on the 
basis of their infrared morphology and color in Spitzer images identifies a signature 
of compact star-forming regions shared across a broad range of luminosities, and by 
inference, masses. We discuss the origin of their striking mid-infrared appearance, and 
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suggest that future studies of the yellowball sample will improve our understanding of 
how massive and intermediate-mass star-forming regions transition from compact to 
more extended bubble-like structures. 

Subject headings: ISM: bubbles—stars: formation—stars: massive—stars: pre-main 
sequence—stars: protostars 


INTRODUCTION 


The Milky Way Project (MWP; ISimpson et al.ll2012l ) is one of a suite of highly productive on¬ 
line citizen sc ience initiatives in the Zooniverse, dev e Mped and maintain ed by the Citizen Science 


Alliance (e.g. iLintott et al.l l2008l : ISmith et ah 


2011 


Fortson et al.l l2012l l . The first implementa¬ 


tion of the MWP utilized archived Syitzer GLIMPSE/MIPSG AL images ( Reniamin et al.l 1200, 'll : 
Mizuno et al.l [20081; iGhurchwell et al.l l2009l : iGarev et al.l [2009|) to study star formation over one 
third of the Galactic plane via the categorization of infrared “bubbles”, which are characteristic 
of H II regions and their associated photodissociation regions (PDRs). The PDRs are prominent 
in the IRAC 8 /rm band, which traces emission from polycyclic aromatic hydrocarbons (PAHs), 
while MIPS images of bubble interiors often show 24 /rm emi ssion, which is likely associated with 
thermal emission from dust grains within H II regions (e.g., IWatson et al.ll2009l i. The principal 
task of citizen science volunteers entailed using an ellipse-drawing tool to mark the sizes, orienta¬ 
tions, ellipticities, and thickness es of H II region/PDR feature s. The first data release expanded 
the previous bubble catalogs of IGhurchwell et al.l (j2006l . l2007li by nearly an order of magnitude 
(jSimpson et al.ll2ni2l ). Furthermore M WP citizen science cla- ssifications have been used as training 
sets for a machine learning algorithm (jBeaumont et al.l 120141 ). 


Demonstrating the serendipitous nature of citizen science efforts, volunteers went beyond their 
assigned tasks and started tagging and discussing, using the MWP ‘Talk’ interface, compact yellow 
objects (“yellowballs”) in the GLIMPSE/MIPSGAL images shortly after the MWP opened to the 
public. The term is descriptive of their color and compact appearance in the GLIMPSE/MIPSGAL 
images, which use a 4.5 /rm (blue), 8 /rm (green), and 24 ^m (red) representative color schem^. 
In total 928 yellowballs were identified by MWP participants (see Table [T]) . Most yellowballs ap¬ 
pear in three types of environments: (1) as isolated objects in filamentary infrared dark clouds 
(IRDGs), along with bright 24 um point sources that are typically associated with embedded mas¬ 
sive protostars fe.g.. iRathborne et al.ll2niril : iBattersbv et al.ll2^14^ : (2) clustered at the intersection 
of bipolar bubbles, which have been associated with outflows from massive protoclusters, and are 
often perpendicular to filamentary IRDGs; and (3) in bubble hierarchies, often along the rims of 
large bubbles, many of which have been associated with known H II regions. 

In Sections [2] and [3] we use the spatial distribution of yellowballs, cross-matches with existing 


^All MWP images viewed by volunteers are available at http://mwp-milkman.herokuapp.com/ 
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catalogs of star-formation tracers, and mid- and far-infrared photometry to show that the yellow- 
balls identified in the MWP are a collection of objects tracing a compact, dense phase of massive 
(O- and B-type) star formation. This includes a mix of compact and ultra-compact H II regions as 
well as analogous regions for less massive B stars. We then examine in Section 0] the origin of the 
striking mid-infrared appearance of yellowballs and discuss how they fit into our current picture of 
massive star formation. Finally, a summary of our findings and our conclusions are presented in 
Section [3 


2. THE YELLOWBALL - STAR-FORMING REGION CONNECTION 


2.1. The Galactic Distribution of Yellowballs 


Figure [T] shows the Galactic longitude and latitude distribution of yellowballs, along with 
the distribution of MWP bubbles, and young Red MS X Source (RMS) objects, which have been 


Lumsden et al. 

2013: 

Urouhart et al. 

2014) 


Each histogram was normalized to its largest bin value and offsets were applied in order to 
facilitate comparison. It is readily apparent that the yellowball distributions have many similarities 
with the MWP bubble and RMS object distributions. In particular, the mean of the distribution in 
Galactic latitude is slightly below the Galactic mid-plane (—0?08), and it also peaks slightly above 
the Galactic mid-plane (in the 0?0 - 0?1 bin). In longitude we see common features associated with 

Galactic structure, such as peaks around I = -|-30° —1-40° and I = —20°-30°, as well as a steady 

decrease in the number of sources beyond I ~ -|-40°. In short, the Galactic longitude and latitude 
distribution of yellowballs supports them being a population of sources associated with Galactic 
star formation. 


2.2. Association of Yellowballs with Dense Clumps/Cores 



to the average size (diameter) of yellowballs (0f4, a ~ 0f3) we find that 245 (49%) of the 502 
yellowballs found within the BGPS survey region have a BGPS association. The ATLASGAL 
survey encompasses the entire MWP region, and, using the same match distance, we find 524 
(56%) of the 928 yellowballs have matches with an ATLASGAL compact source. The catalog 
identifications for associated BGPS and ATLASGAL sources are listed in Columns 2 and 3 of 
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Table [2j 

For comparison, cross-matches done between 20 randomly generated datasets (with statistically 
the same latitude and longitude distribution as the yellowball sample) had only 6 ± 2 (1 ± 0.4%) 
and 10 ± 2 (1 ± 0.2%) BGPS and ATLASGAL associations on average respectively. We conclude 
that any contamination from random associations is minimal, and the proximity of yellowballs to 
regions of dense molecular gas is again consistent with what would be expected for a population of 
objects associated with star formation activity. 


2.3. Association of Yellowballs with H II Regions 


We performed a sim ilar analysis to cross-m atch the yellowball positions with the WISE catalog 
of Galactic H II regions (jAnderson et al.ll2014l i. We found that 599 (65%) yellowballs have matches 
to within 0l4 (see Column 5 of Table [2]). Entries in the WISE H II region catalog were selected 
by employing mi d -infryed cr iteria similar to those used to identify “bubbles” in the MWP and 
Churchwell et al.l ([200g, 120071 ) catalogs hence it is likely not complete for very compact sources. 
Their resulting catalog differentiates between four classes of objects (K, G, C, Q): known H II 
regions, which have measured Radio Recombination Line (RRL) or Ha spectroscopic emission 
(K); grouped H II regions, where candidates are associated via positional correlation with known 
multiple H II regions (G); candidate H II regions, which have characteristic H II region mid- 
infrared morphology spatially coincident with detected radio continuum emission but lack RRL or 
Ha observations (G; follow-up observations by Anderson and colleagues suggest that essentially all 
objects in group G are bona fide H II regions); and radio quiet objects, which may contain only 
intermediate-mass stars, or be H II regions in either early or late stages of evolution (Q). However, 
since most of the radio quiet objects have small angular sizes and correlate with c old dust, they 
are probably in the earliest phases of H II region evolution (jAnderson et al.M2ni4l ). Yellowballs 
with WISE H II region catalog matches span all four categories, with associations as follows: 185 
K (31%), 62 G (10%), 144 G (24%), and 208 Q (35%). The WISE source class for each matched 
yellowball is shown in Column 6 of Table [2j 


2.4. Association of Yellowballs with Red MSX Sources 


The RMS catal og is the largest stati stically selected catalog of young massive protostars and 
H II regions to date (jLumsden et al.ll2013l ). It correlates spectral information across a wide range of 
wavelengths from the near-infrared to the radio regimes, uses rigorous color criteria to classify dif¬ 
ferent types of objects, and, where possible, includes kinematic distance estimates and bolometric 
luminosities for catalog entries. The final catalog lists II categories of objects, including five cate¬ 
gories associated with evolved star groups (generic evolved stars, planetary nebulae (PNe), proto¬ 
planetary nebulae, OH/IR stars, and carbon stars); four categories of young objects (YSOs, H II 
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regions, H II/YSO, and diffuse H II regions); and two categories of ambiguous objects (young/old 
sources and other). The catalog is thought to be complete for the detection of a BO V star at the 
distance of the Galactic center, although inclusion in the catalog is dependent upon detection of 
the source by MSX and specific RMS color criteria (such as rising flux toward longer wavelengths 
in the MSX bands). 


Lumsden et al.l (120131 1 estimate that 95% of the Galactic ultra-compact H II regions were 


detected, but more than 50% of the larger compact H II regions may be missing, although some of 
these might have been classified as “diffuse H II regions” in the RMS catalog. Therefore, the RMS 
classification “H II region” should generally be interpreted as “ultra-compact H II region”, while 
“diffuse H II region” may in fact include compact H II regions. 


To investigate association of yellowballs with RMS objects, we first cross-matched the yel- 
lowball list with the RMS catalog using a cross-match distance corresponding to the size of each 
individual yellowball (Urquhart, private communication). Of the 825 yellowballs that overlap the 
RMS catalog (103 yellowballs lie within 10° of the Galactic center, and would therefore not be 
included in the RMS catalog), 282 (~ 34%) have RMS matches (see column 4 of Table [2]). Of 
the 282 yellowballs with RMS matches, 155 (19%) are positionally coincident with RMS objects 
to within 5", indicating the yellowball is the main contributor to the luminosity of these sources. 
Of these 155 yellowballs, 18 are classified as “Diffuse HII region”, 4 as “HII/YSO”, 17 as “YSO”, 
and the remaining 116 as “HII region” in the RMS catalog. No yellowballs are associated with 
any of the five evolved star groups or ambiguous objects in the RMS catalog. All but 17 of the 
155 yellowballs have bolometric luminosity estimates; these span a broad range from 3.30x10^ - 
1.8x106 Lo (see Figure [2]), with the vast majority of the yellowballs having luminosities expected 
for regions of massive star formation. 


Figure [3] shows the distribution of physical sizes for the same sample of 138 yellowballs. For 
each yellowball an effective angular size was calculated using the average of Al and A6 from Tabled] 
then converted to a physical size using distances from the RMS catalog. It should be kept in mind 
that these sizes are strict upper limits to the sizes of any ionized regions within these objects for two 
reasons: (1) the size of each yellowball corresponds to the size of the user-drawn rectangle enclosing 
the yellowball, and (2) the yellowball is tracing the maximum extent of the PDR associated with 
each source so any enclosed H II region is necessarily smaller (see Section dj) . This suggests that 
many of the smaller yellowballs are likely to contain ultra-compact H II regions or perhaps even 
hypercompact H II regions, ident ified with the earliest manifestations of ionized gas around young 


massive stars (|Franco et al 


200 ' 


3). 


Figure d) presents two images from the MWP containing yellowballs with and without RMS 
counterparts (coincident to within 5''). The top panel shows the central portion of W 33, with 
a cluster of yellowballs located between two large infrared bubbles. Yellowballs 541 and 880 are 
classified as “YSO” and “diffuse H II region”, respectively, in the RMS catalog, while yellowballs 49 
and 542 are classified as H II regions. The latter are unresolved near the center of the image, with 
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a “lima-bean” morpholo gy. Yellowball 450 is a known H II region, first detected as a 20 /rm s ource 
by iDvck &: SimonI (jl977l i and characterized at cm radio wavelengths by iHaschick &: Hoi ([l983|) , but 
it was not included in the RMS catalog because it was not included in the MSX catalog. 


In the lower panel, yellowballs 78 and 643 both have RMS counterparts, and are classified as 
an “H II region” and “YSO” respectively. In contrast, the other yellowballs in the field (63 and 
697) fail the color criteria for inclusion in the RMS catalog: both have MSX F(8 //m) > F(14 
/im), which is likely due to them having a large PAH ionization fraction (see Section [4] for more 
discussion). 

The MWP yellowball catalog provides a useful compliment to the RMS survey. It similarly 
highlights regions of massive star formation, but our results suggest that the selection of objects 
purely on the basis of their infrared morphology and color in the higher-resolution Spitzer images 
(cf. MSX) identifies a signature of compact star-forming regions shared across a broad range of 
luminosities, and by inference, masses. 


3. PHOTOMETRY 


We performed an infrared color analysis in order to explore further the association of yellow¬ 
balls with star-formation activity apparent both by visual inspection of MWP images and by the 
numerous catalog cross-matches outlined in the previous section. To rule out the possibil i ty tha t 
these compac t objects might be PN e. we used color criteria developed bv lAnderson et al.l (j2ni2l ). 
In particular. lAnderson et al.l (120121 ) showed that the 12 fiui to 8 /um flux ratio, derived from WISE 
and Spitzer IRAC data respectively, is a robust discriminator for separating H II regions and plan¬ 
etary nebulae (PNe). The IRAC 8 band is very sensitive to emission associated with the PAH 
bands centered at 7.7 and 8.6 /xm. In contrast, although the WISE 12 fim band does cover the 7.7, 
8.6, 11.3 and 12.7 /xm PAH bands, its normalized response is relatively low (~ 0.4 — 0.8) at these 
wavelengths, and the peak sensitivity of the band lies at 14.5 /xm, longwards of most of the PAH 
emission. Given this combination of PAH and filter properties, H II regions tend to have a lower 
12/8 flux ratio than PNe because their mid-infrared emission is dominated by stro ng, broad PAH 


featu res, whereas PAH features in PNe tend to be narrower and relatively weaker (iBregman et al 
1989l j. 


3.1. Selection Criteria & Technique 

We performed aperture photometry on a representative sample of 183 yellowballs using Spitzer 
8 /xm &: 24 /xm and WISE 12 /xm band data. This sample was chosen using all yellowballs coincident 
with RMS sources (to within 5") that were not saturated in MIPSGAL images at 24 /xm (N=81), 
and a comparable number of yellowballs without RMS matches (N=102). The latter sample was 
chosen by selecting yellowballs with the highest (N=52) and lowest (N=50) hit rates, where the 
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hit rate is defined by the ratio of the number of times a yellowball was identified by a volunteer 
to the number of times it was viewed. The rationale for choosing this sample was to allow us to 
investigate any potential differences based on hit rate. Only objects viewed > 50 times, which have 
hit rates >0.1, are included in the list of 928 yellowballs. 


Aperture photometry was done using the IDL-based imview program (jHiggs et al.lll997l ). For 
extended infrared sources at low Galactic latitude, especially those associated with star-formation 
activity, background estimation and subtraction is the most important sou rce of photometric erro r 
as the infrared emission in the Galactic plane is highly structured (e.g., iFich &: TerebevI Il996l f . 
IMVIEW allows the user to select points that define the shape of the photometric aperture, and that 
constrain the surface fit used to estimate the background. The background fit can then be varied 
by selecting different background points and/or using different interpolating functions. For each 
source the average flux, using two different apertures and two different surface fits, was obtained 
and the standard deviation of the measurements was used as a measure of the uncertainty. 


3.2. Results 


The average flux density and uncertainty at 8, 12, and 24 /rm for each yellowball is reported 
in columns 2 ~ 7 of Table [3j The table is divided into three subsections corresponding to the 
RMS-match sample, and the high and low hit rate non-RMS samples. In addition, Golumns 8-11 
of Table [3] indicate whether the source has a cross-match in another catalog, and gives the WISE 
catalog class if applicable. 


In Figure[5]we show a histogram of the resulting log(Ti 2 /F 8 ). The average valu e for this sample 


(—0.1 9), and the average color for the sample of WISE H II regions (—0.09) from lAnderson et al 


(j2012l ) are indicated. The photometric uncertainties are typically about half of a bin width (~ 0.05), 
and the entire sample (except for one yellowball with a morphous bound a ries) satisfies the robust 
F(12 /im)/F(8 /rm) flux ratio criterion determined bv lAnderson et al.l (120121 ) to separate H II 
regions from PNe at -1-0.3. We note that the RMS and non-RMS subsamples do not have identical 
distributions, specifically the average color of the non-RMS yellowballs (—0.23) is more negative 
than the average of the RMS-matched subsample (—0.13), and a KS-test shows that the two samples 
are significantly different (significance level p = 0.001). We also find that, within the non-RMS 
sample, the high hit rate sample has an average color (—0.28) that is more negative than the low hit 
rate sample (—0.17). Again, a KS-test shows the two samples are signficantly different (significance 
level p = 0.001). We discuss the likely cause of these color differences in the next section. 


We found that yellowballs have on average F 24 ~ 3.5^8 {a = 2.5, median F 9 A 2.9^8). 
is consistent with the average H II region spectrum shown in lAnderson et al.l (120121 ) . 


This 
















4. DISCUSSION 


The distinct mid-infrared appearance of the yellowballs in the GLIMPSE/MIPSGAL images 
used in the MWP is not primarily due to the rough equality of the 8 /um and 24 |Um fluxes 
mentioned in Section [3l as this is expected for all H II regions, but comes about because the 
emission is spatially coincident. This spatial coincidence is expected in the early stages of the 
evolution of H II regions/PDRs. Models predict that large PDRs will form around any initial 
(dust-filled) ionized region, and that the maximum size of the PDR, and the time at which the 
maxi mum size is obtained, are b oth relatively insens i tive to stellar luminosity (iRoger fc Dewdnev 
19921 ) ■ For example. Figure 9 of lRoger fc DewdnevI (Il992l ) shows that the PDR associated with 
an H II region expanding into a molecular cloud having densities between 300 to 3000 cm“^ will 
reach a maximum thickness of order 0.1 to 1 pc on timescales of 10 ^ — 10 ^ years for a wide range 
of stellar luminosities. This size range of the model PDR, and the insensitivity of the model to 
stellar luminosity, are both consistent with the derived yellowball sizes and the range of luminosities 
presented in Section 2.4. 

As the region evolves the ionized region is expected to catch-up to the photo-dissociation front 
resulting in a thin, shocked H I region/PDR surrounding the H II region. Much of the dust will be 
remove d from the ce ntral portion of the H II region via the action of radi ation pressure and stellar 
winds (jDrainell201lh . and PAHs will be destroyed within the ionized gas (jGiard et al.lll994l '). This 
leads to a clear spatial separation between the Tg emitting region (PAH-rich, PDR) shown as green 
in the MWP images, and the T 24 emitting region (depleted interior hot dust, perhaps resupplied by 
the erosion of denser clumps in the region; lEverett &: Churchwellll2010l f shown as red in the MWP 
images. 

We noted in Section [3] that the yellowball sample has an average log(Fi 2 /T 8 ) color that is 
more negative than the average color for the general H II region sample. This is most likely 
due to the fact the yellowball sample contains a higher fraction of compact objec t s with a large 
PAH ionization fraction (Ybarra 2014, private communication). iRoelfsema et al.l (jl99m showed 
that compact H II regions had a much higher 7.6/11.2 PAH intensity ratio (ranging from 5 - 
11) compared to more evolved H II regions (~ 3). This effect is related to the PAH ionization 
fraction because the strength of the PAH bands around 7-8 /um is highly sensitive to the ionization 
state of the PAHs (becoming much stronger in ionized PAHs), where a-s the strength of P AH bands 
at longer wavelengths are minimally affected by the ionization state (|Draine fc Lil 120071 ). Thus we 
would expect the log(Fi2/T8) color to become increasingly negative with increasing PAH ionization 
fraction. 


This may also be the cause of the shift in the average log(Fi 2 /T 8 ) color from —0.13 for RMS- 
matched yellowballs to —0.23 for yellowballs with no RMS matches. The median size of the non- 
RMS matched yellowballs was 0l29 compared with 0l46 for the RMS-matched sources. Assuming 
that the non-RMS and RMS-matched yellowballs have a similar distribution in distance, this implies 
that the non-RMS sources are more compact on average. 
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Similarly we observed a difference in the average color of high (—0.28) and low hit rate (—0.17) 
yellowballs without RMS matches. The median size of the high hit rate sample is 0'.27 compared 
with 0f35 for the low hit rate sample. For the entire yellowball sample the correlation between hit 
rate and angular size is very weak (correlation coefficient r = —0.08). This is not surprising, as 
there are clearly other factors such as the proximity of the yellowball to other interesting objects, 
and the overall complexity of the held, that inhuence hit rate. Given this weak correlation, if the 
two samples have a similar distribution in distances, then the high hit rate sample will have a higher 
proportion of physically compact objects leading to a more negative average log(Ti 2 /F 8 ) color. 


As touched on in Section [2l it is clear that the RMS catalog does not include some yellowballs 
that are dehnitely star-forming regions. To explore this further we cross-matched the 646 non-RMS 
yellowballs with the MSX catalog using a search radius of lO'k Of these 199 (31%) did not have an 
MSX catalog entry, 255 (39 %) had poor quality (S/N < 5) in band E (21.3 um), and the remainin g 
192 (30 %) failed one or more of the MSX color cuts described in detail in iLumsden et al.l (j2013ll . 
This last subset is particularly interesting as it probably includes many compact objects with strong 
PAH emission. 


The catalog cross-matches presented in this paper offer opportunities for a multitude of pos¬ 
sible explorations. For example, consider the 138 yellowball-RMS matches that have distance and 
luminosity estimates. Of these objects, 35 have BGPS counterparts and 114 have ATLASGAL 
counterparts, indicating that nearly all yellowballs with RMS entries are associated with dense 
gas. Objects in this group span the range of luminosities seen in Figure 2. On the other hand, 
of the 208 yellowb a lls th at are classified as radio-quiet objects in the WISE H 11 region catalog 
of Anderson et al. ( 2014 1. only 96 (46%) have BGPS or ATLASGAL counterparts and only 25 
(12%) have RMS associations. Twenty-two of the 25 RMS-associated radio-quiet yellowballs have 
luminosity estimates, all of which are < 5 x IO^Lq (< BO ZAMS equivalent), suggesting that many 
low-luminosity yellowballs were missed by these surveys. 

The most straightforward interpretation of these low-luminosity yellowballs is that they are 
associated with the formation of mid- to late- B stars, which are expected to have large, relativel y 


long-lasting PDRs, combined with small, weak H II regions (jKerton 


2002 


Lundquist et al.l 120141 1. 


The long lifetime of the PDRs in this case could explain why the percentage of yellowballs with 
dense gas associations is lower than that found for the RMS-matched sample. Intriguingly though, 
it is possible that some of these sources are massive protostars in a pre-ultra-compact H II region 
stage; such protostars can exhibit low overall luminosities due to low core temperatures while 
at t he same time having mid-infrared emission t hat would be visible as a yellowball (e.g., Mol 
160, iMolinari et al.l 120081 : IWolf-Ghase et al.l l2012l ). Further high-resolution and high-sensitivity 
observations at infrared and radio wavelengths would help to distinguish between these two options 
by, for example, detecting outflows expected for massive protostars, or by detecting the small H II 
regions expected to be associated with the B-type stars. 
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5. SUMMARY & CONCLUSIONS 


We have presented multiple lines of evidence that yellowballs are a mix of compact star-forming 
regions, including ultra-compact and compact H II regions, as well as analogous regions forming 
less massive stars. Visual inspection of MWP images indicates these objects are typically found 
in IRDCs and/or in bubble hierarchies, and their distribution in Galactic longitude and latitude 
mirrors the distribution of MWP bubbles and young RMS objects. Color analysis of yellowballs 
using Spitzer and WISE data indicates that yellowballs occupy regions of infrared color space that 
include H II regions and exclude evolved objects such as PNe. Cross-matching the yellowballs 
with the ATLASGAL, BGPS, WISE-H II region, and RMS catalogs, indicates that the majority 
of these objects are unambiguously associated with dense molecular clumps and other signposts of 
star formation. No yellowballs are associated with any of the evolved star categories in the RMS 
catalog. 

We examined the luminosity and physical size distribution of a sample of 138 yellowballs 
with < 5" positional associations with RMS catalog sources. These objects span at least 3 orders 
of magnitude in luminosity, and have luminosities consistent with massive star-forming regions. 
Typical sizes for these yellowballs are comparable to compact H II regions (sub-parsec in size); 
however, we expect that these sizes are upper limits on the extent of the ionized gas as the size 
includes emission from the surrounding PDR. The great majority of these sources are ultra-compact 
H II regions or even younger/denser objects. 

The origin of the distinct yellow color of these objects in the GLIMPSE/MIPSGAL MWP im¬ 
ages is the cospatial emission from PAHs and dust. This is expected for the earliest stages of massive 
star formation when PDRs will be at their thickest extent and various dust-clearing/destruction 
mechanisms within the ionized gas have had minimal time to act on the dust distribution. 


Yellowballs are analogous to the “green peas” of the archetype Galaxy Zoo project, in the 
sense that they represe nt a class of objects identified and recorded by citizen science volunteers 
bnod b Just as the discovery a nd subsequent studies of green peas ha.ye yie lded 


(Gardamone et al, 


critical new insights into the evolution of galaxies (lAmorin et al.ll20I2l : iGhakrakorti et al.ll20I2l ). we 
expect that future studies of the yellowball sample will improve our understanding of how massive 
star-forming regions transition from compact embedded stages (e.g., massive protostars and ultra¬ 
compact H II regions) to more evolved H II regions. Similarly, those yellowballs associated with 
slightly lower-mass star formation will provide us with a comparable view of how these objects 
transition from highly embedded stars/clusters to larger bubble structures. 
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that a high PAH ionization fraction could produce more negative log(Ti 2 /T 8 ) values based on his 
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James Urquhart for his assistance with cross-matching the yellowball list and RMS catalog, and 
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Fig. 1.— Galactic longitude (upper) and latitude (lower) distribution of yellowballs (YB). The 
distributions are similar to the other tracers of Galactic star formation shown; young RMS objects 
(RMS), MWP large bubbles (MWP-L), and MWP small bubbles (MWP-S). Note RMS does not 
cover |/| < 10°. Each histogram has been normalized to its largest bin value and an offset of 0, 1, 
2, and 3 has been added to the MWP-L, MWP-S, RMS and YB histograms respectively. 
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Fig. 2.— Yellowball Luminosity. Luminosities from iLumsden et al.l (|2ni,'ll ) are shown for a subset of 
138 yellowball s with RMS asso ciations (coincident to within 5''). For reference the ZAMS luminosity 
calibration of ICrowtherl (|2005l i is shown along the upper axis. 
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Fig. 3.— Yellowball Angular Size and Distance. Points are plotted for the 138 yellowballs (YB) 
coincident to within b" with RMS sources. The horizontal dashed lines show the and 95**^ 
percentiles for the entire yellowball sample (e.g., 95% of all yellowballs have an angular size < ih 
The diagonal solid lines represent objects with a constant physical size (as labeled). For reference 
the nominal MIPSGAL and GLIMPSE resolutions are shown as horizontal dash-dot lines. 
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Fig. 4.— MWP GLIMPSE/MIPSGAL images, 4.5 fivu (blue), 8 /xm (green), and 24 /xm (red), 
showing yellowballs with (in white) and without (in yellow) RMS counterparts. The top panel 
shows part of W 33, with a cluster of yellowballs located between two large bubbles. The lower 
panel shows two isolated yellowballs and two yellowballs associated with a larger H II region (see 
the main text for details). The Galactic coordinates of the image center and the angular size of the 
image are shown at the top of each frame. 
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Fig. 5.— WISE and Spitzer photometry for a sample of 183 yellowballs (see the main text for 
selection criteria). The dotted line shows the avera ge color (—0.19), which can be compared with 
the average color of the sample of H 11 regions from lAnderson et al.l (120121') shown as a short-d ash 
line (—0.09). The long-dash line shows the H 11 region - PNe cutoff from lAnderson et al.l (j2012l ') at 
-1-0.3. The grayed region shows the distribution of yellowballs without RMS associations (N=102; 
average color = —0.23). 
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Table 1. Yellowballs Identified by MWP Participants 


Short 

ID 

MWP 

ID 

1 

(deg) 

b 

(deg) 

AP 

(arcmin) 

Ab^ 

(arcmin) 

Hit 

rate*’ 

1 

MWP1G311565-F00230Y 

311.565 

0.230 

0.160 

0.150 

0.78 

2 

MWP 1 G340572-h00360Y 

340.572 

0.360 

0.349 

0.338 

0.78 

3 

MWPlG342062-h00422Y 

342.062 

0.422 

0.496 

0.543 

0.67 

4 

MWPlG032122-h00091Y 

32.122 

0.091 

0.503 

0.500 

0.65 

5 

MWPlG327901-h00154Y 

327.901 

0.154 

0.259 

0.263 

0.65 


Note. — Table [T] is published in its entirety in the electronic edition. A portion is 
shown here for guidance regarding its form and content. 

'^Average size of user-drawn rectangle enclosing the yellowball. 

'^The ratio of number of times the yellowball was identihed to number of times the 
yellowball was viewed. 






Table 2. Yellowball Catalog Matches and MIPSGAL 24 //m Saturation 


Short ID 

BGPSv2'‘ 

ATLASGAL 

RMS 

WISE 

WISE Glass*’ 

MIPSGAL Saturated" 

1 


G311.5638-b0.2298 

N 

N 


N 

2 


G340.5742-b0.3608 

N 

G340.573-b00.359 

Q 

N 

3 


G342.0579-b0.4211 

G342.0610+00.4200 

G342.062-b00.417 

K 

Y 

4 

G032.119-b00.092 

G032.1173-b0.0909 

N 

G032.123-b00.086 

G 

N 

5 


G327.9065-b0.1573 

G327.9018+00.1538 

N 


Y 


Note. — Table[2]is published in its entirety in the electronic edition. A portion is shown here for guidance regarding its form and content. 

Note. — If a catalog cross-match was made the catalog identification is shown. N indicates no cross-match was found. 

'‘Sources not in the BGPS survey area have • • • entries. 

^See text for details. If no WISE cross-match was found • • ■ is shown. 

"Indicates if any portion of the source was saturated (Y) or not (N) in MIPSGAL 24 fim images. 


K) 

O 






Table 3. Yellowball Photometry 


Short ID Fs 

(TFg 

Fi2 

0'Fi2 

F 24 

(7^24 BGPS^ ATLASGAL"- WISE"- WISE Glass'^ 

(Jy) 

(Jy) 

(Jy) 

(Jy) 

(Jy) 

(Jy) 


RMS Matches 


46 

65 

68 

0.44 

2.61 

0.64 

0.05 

0.04 

0.09 

0.31 

2.00 

0.67 

0.05 

0.05 

0.12 

3.23 

6.73 

3.92 

0.22 

0.17 Y 

0.71 

Y 

Y 

N 

Y 

Y 

Y 

G 

G 

G 

No RMS Match - High Hit Rate 

1 

0.05 

0.01 

0.03 

0.01 

0.81 

0.02 

Y 

N 


2 

1.32 

0.05 

0.83 

0.08 

1.24 

0.04 

Y 

Y 

Q 

4 

5.35 

0.22 

3.75 

0.24 

15.04 

0.53 Y 

Y 

Y 

G 

No RMS Match - Low Hit Rate 

847 

0.023 

0.002 

0.011 

0.003 

0.011 

0.002 

N 

N 


848 

0.074 

0.005 

0.07 

0.01 

0.29 

0.03 

N 

N 


850 

0.82 

0.05 

0.67 

0.04 

1.76 

0.07 

N 

Y 

Q 


Note. — Table E] is published in its entirety in the electronic edition. A portion is shown here for guidance 
regarding its form and content. 

‘^Y^Catalog match, N= No catalog match, • • ■ = source not in survey area. See Table [2] for match source 

ID 

’’ • • • is shown for sources without a WISE match. 












